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'sj" ■ Abstract 

o : 

O . It is useful to study the space of all cosmological models from a dynamical systems 

I perspective, that is, by formulating the Einstein field equations as a dynamical sys- 
tem using appropriately normalized variables. We will discuss various aspects of this 

. work, the choices of normalization factor, multiple representations of models, the past 

^ I attractor, nonlinear dynamics in close-to-Friedmann-Lemaitre models, Weyl curvature 

fS| ' dominance, and numerical simulations. 

1 Introduction 

00 : 

^ I The hot big-bang model of modern cosmology is based on the assumption that the large 
O ' scale geometry and dynamics of the universe can be described by an exact FLRW model, 
^ ' with more detailed and smaller scale physical phenomena, such as density fluctuations or 
■ gravitational waves, being described by perturbed FL models, i.e. by solutions of the EFE 
qh| linearized about an exact FL model. 

^ [ We begin with the premise that it is important to study cosmological models more general 

^1 than perturbed FL models. In particular one is interested in the set A^obs, namely, the set 
' of all universes that are compatible with current observations (Ellis 2004). A member of this 
^ . set must have an epoch, possibly finite, during which the model is close to FL; at early or 
later times it may deviate significantly from FL. 

More generally we are interested in the set Ai of all cosmological models, and in the 
relationship between A4, A^obs and A^fl, the set of all FL models. For example, we want 
to study the relation between linearly perturbed FL models and models that are close to 
FL in some well-defined, gauge-invariant sense, and satisfy the exact EFE. Are nonlinear 
effects perhaps significant in close-to-FL models? Another question that takes one outside 
the realm of linearly perturbed FL models is the detailed nature of the generic cosmological 
singularity. 

With the preceding discussion as motivation we consider the class Ai of all cosmological 
models, whose state at time t can be represented by a vector X in a state space S. The 
evolution of a model universe is then described by a curve X = X(t) in S, called an orbit, that 
will be a solution of a system of first order autonomous evolution equations and constraints 
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of the form 



c(x,a,x,---) = o, 

where di denotes partial differentiation with respect to the spatial coordinates x\ and ■ ■ • 
denotes possible higher order spatial derivatives. The goal is to choose state variables that 
remain bounded during the evolution of the model and a time variable t such that t — > — oo 
at the initial singularity and t — ^ +00 at late times. 

Within this framework, classes of models with symmetries will be represented by invariant 
subsets of S, the most important being the subset iSfl that describes the FL cosmologies. 
One expects that the orbits of close-to-FL models will shadow orbits in iSfl- 

The appropriate mathematical vehicle for implementing the above program is the or- 
thonormal frame formalism^, since it expresses the EFE directly as first order (in time) 
autonomous evolution equations. It is also necessary to introduce a process of normalization 
in order to create bounded variables. A choice that has proven effective is the so-called Hub- 
ble normalization. This approach to the study of cosmological dynamics has been discussed 
in some detail in Wainwright & Ellis 1997^, and more recently by Coley 2003, with emphasis 
on spatially homogeneous models. Since then the approach has been extended to models 
without symmetry by Uggla et al. 2003. 

The plan of this paper is as follows. In Section |21 we describe the principal features 
of the Hubble-normalized state space, and in Section El we discuss close-to-FL models and 
the related notion of isotropization. In Section |3] we give an overview of the dynamics of a 
special class of models in the cosmological hierarchy, the nontilted spatially homogeneous 
cosmologies. In Section El we describe a candidate for the past attractor in the Hubble- 
normalized state space, and mention some recent numerical simulations which support this 
proposal. In Section IHl we briefly introduce a modification of Hubble normalization which 
provides a description of nontilted spatially homogeneous models that undergo recoUapse. 
Section [71 contains our concluding remarks. 

2 The Hubble-normalized state space 

In this section we give the motivation for introducing Hubble-normalized variables, and 
describe the principal features of the resulting Hubble-normalized state space and evolution 
equations. 

^This formalism was introduced in relativistic cosmofogy by Ellis 1967, among others, and an extended 
version was given by van Elst & Uggla 1997. 

^We shall abbreviate this reference to WE in this paper. 
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2.1 Hubble-normalized variables 



One of the rationals for using dynamical systems methods in cosmology is the hope that one 
can describe the evolution of cosmological models near the initial singularity by means of 
a past attractor of a cosmological dynamical system, and the dynamics at late times in a 
particular epoch by means of a future attractor. 

In order to formulate the Einstein field equations as a dynamical system it is clear that one 
has to normalize the variables, since near the initial singularity physical variables typically 
diverge and at late times typically tend to zero. Physical considerations suggest that in a 
cosmological setting, normalization with the Hubble scalar of the fundamental congruence is 
an appropriate choice. Firstly, consider the density parameter Qm, which plays a fundamental 
role in cosmology, in that it measures the influence of the matter on the dynamics. In 
geometrized units, it is defined to be the matter density p divided by 3i7^, where H is the 
Hubble scalar. Secondly, the extent to which the overall expansion of the universe is close to 
isotropy is measured by the ratio of the shear (the trace-free part of the expansion tensor) 
to the Hubble scalar. 

The importance of Hubble normalization appears to have first been mentioned in the 
literature by Kristian & Sachs 1966 (see page 398) in connection with their general analysis 
of the potential constraints that observations of distant galaxies can place on the geometry 
of spacetime. The role of Hubble-normalized variables emerges again when one analyzes 
the potential constraints that arise from observations of the CMB (Maartens et al. 1995, 
1996). In both analyses one finds that bounds are placed on Hubble-normalized physical 
and geometrical quantities. Our principal interest, however, lies in using Hubble-normalized 
variables as dynamical variables, which dates back to Collins 1971, who gave a qualitative 
analysis of the dynamics of some special classes of cosmological models. 

Working within the orthonormal frame formalism we introduce Hubble-normalized vari- 
ables and write the EFE as first order evolution equations in the Hubble-normalized state 
space. For simple classes of ever-expanding models, for example open FL models, and spa- 
tially homogeneous models of Bianchi type I, the Hubble-normalized state space is bounded. 
Near the singularity physical variables such as the matter density diverge, but the Hubble 
scalar also diverges and at such a rate that the corresponding Hubble-normalized quantity, 
the density parameter Qm = f,„ , remains bounded. In contrast, if H is close to zero, 

then the other physical variables are equally close to zero, so that the Hubble-normalized 
quantities are bounded. 

For more general models, however, it turns out that this simple picture is not valid 
- the Hubble-normalized state space is unbounded. The primary reason is that the Hubble- 
normalized Weyl curvature tensor can assume arbitrarily large values. All available evidence, 
however, suggests that the Hubble-normalized variables are bounded into the past, i.e. on 
approach to the initial singularity. In addition, if there is a positive cosmological constant 
there is strong evidence that the Hubble-normalized variables are also bounded into the 
future. Thus, even though the Hubble-normalized state space is unbounded, one expects 
that the evolution equations will admit a past attractor and a future attractor. 
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2.2 Asymptotic regimes 

When one formulates the EFE as a dynamical system, one uses a time variable t that 
potentially assumes all real values. The asymptotic regimes are then defined by the limits 
t — > — oo and t — > +00. Cosmologists model the physical universe as a sequence of epochs 
in time, an epoch being identified by which source term is dynamically dominant. A typical 

succession of epochs is 

i) inflationary, dominated by a scalar field (p, 

ii) radiation-dominated, p = \p-, 

iii) matter-dominated, p = 0, 

iv) accelerating, dominated by a cosmological constant A > 0. 
With each source term is associated a Hubble-normalized energy density, 

^rnt 

and a particular epoch is defined as the time interval during which a particular Q is dominant. 

In analytic work one typically assumes the presence of one or two source terms^ and 
one investigates the asymptotic regimes, which are described mathematically by the past 
attractor and future attractor of the dynamical system. Thus, for example, the past asymp- 
totic regime of the radiation-dominated epoch would physically coincide with the end of the 
infiationary epoch and the future asymptotic regime would coincide with the beginning of 
the matter-dominated epoch. 

2.3 The gravitational and matter variables 

A cosmological model is a spacetime with a preferred timehke vector field u whose metric is 

a solution of the EFE with appropriate matter/energy content. It is assumed that there is 
an epoch during which the model is expanding, i.e. the preferred vector field satisfies 

H := IVaU" > 0, 

where H is referred to as the Hubble scalar. We assume furthermore that u is the normal 
vector field to a family of spacelike hypersurfaces t = constant, thereby defining a time 
variable. 

We introduce an orthonormal frame {6^}, with eo = u, and choose the spatial frame 
vectors to be Fermi-propagated. Within this framework, the gravitational field variables 
are the commutation functions ■j'^ab associated with the frame, defined by 

cosmological constant can be incorporated with little increase in complexity in the evolution equations. 
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They may conveniently be expressed in terms of geometric quantities 



^a/s}, (1) 

according to 

7°o/3 = - Hdfs", 7%„ = Ua, Yap = e^fiu n'"' + 5^ - ap 5,,^, 

(see for example WE, page 32). Here iia is the acceleration and is the shear of the 
congruence u, while Cq, and n^p determine the curvature of the spacelike hypersurfaces 
(WE, pages 18, 19, 34). 

As regards matter/energy content we consider a cosmological constant A and a perfect 
fluid with 4-velocity vector field u and linear barotropic equation of state p = (7 — l)p, with 
1 < 7 < 2. In the orthonormal frame formalism the perfect fluid is described by its energy 
density p relative to eg and its velocity relative to eg, i.e. the projection of u orthogonal to 
eo (see Uggla et al. 2003, equations (2.1)-(2.6) for details). So the gravitational variables ((H) 
are augmented by the matter variables 

{p, A}. (2) 
The velocity is dimensionless and satisfies 

< v^v"" < 1. 

We now introduce Hubble-normalized variables^ and differential operators according to 

a,, n^p}/H, {fi^, n^} := {p, A}/(3if2), (3) 

It is also necessary to introduce the derivatives of the normalization factor, namely the 
deceleration parameter q and the spatial Hubble gradient Va, defined by 

q + l:= -^doH, r„ := -^d^H. (4) 

At this stage we introduce the separable volume gauge (Uggla et al. 2003), which is charac- 
terized by 

and the existence of local coordinates t, such that 

do = dt, da = Eadi , 



^Hubble-normalized variables have also been defined for cosmological models with source terms other 
than a perfect fluid with linear equation of state and a cosmological constant, e.g. magnetic fields (LcBlanc 
1998), scalar fields (Coley & van den Hoogen 2000) and brane matter fields (Coley etal 2004). 
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with the preferred hypersurfaces being given by t = constant. Here dt and di denote partial 
differentiation with respect to t and x*. The deceleration parameter q can then be expressed 
in terms of the variables in Q and the differential operator d^, using the Raychaudhuri 
equation. 

With the above choices of gauge the Hubble-normalized state vector becomes 

X = [Ea , Tq., Aa, Nap, ^m, ^A; ^a) ■ 

The EFE, the Jacobi identities and the commutators lead to a system of evolution equations 
and a system of constraints for the components of X. In order to describe the structure of 
these equations, it is convenient to decompose the state vector X as follows: 

X={Ea\raf(SY, (5) 

where 

Y = (Sa/3, Aa, Na/3, ^m, ^A, Va) ■ (6) 

The evolution equations have the following form: 

ai£;,* = (g5/-S/)V (7) 
atr„ = (g5/-S/)r^ + E„^9,g (8) 
dtYA = Fa{Y) + F^^"(Y)E„*9,Fs + FA''^{Y)Ea'd,rp + F^"(Y) r« , (9) 

where the deceleration parameter q is given by 

g = + i [(37 - 2) + (2 - 7H ^^ _ 
^ 1 + (7 — l)t>^ ^ 

There is also a set of constraints that can be written symbolically as 

= C(X, 9,X). (11) 

One of the constraints, the so-called Gauss constraint, is of particular importance, and we 
thus give it specifically^: 

+ f^A + ^^fc + = 1, (12) 
where Y? is the Hubble-normalized shear scalar given by 

Y? = (13) 

and flk is the Hubble-normalized spatial curvature® of the hypersurfaces t = constant, given 



^Lim et al. 2004, equation (2.27). 

^rife is defined by Qk — ~q§^i where is the curvature scalar of the metric induced on the hypersurfaces 
t — constant. 
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by 

= -K^a - r^)A" + \N^pN^^ - U^an' + A^A". (14) 

The detailed form of equations (f7j)- (fTT|l is given in Lim et al. 2004 (equations (2.18)- 
(2.30) with i?" = 0). What concerns us here is the overall structure of this system. We make 
a number of observations. 

i) The equations are first order evolution equations in time, but do contain second order 
spatial derivatives, namely the second derivatives of the spatial Hubble gradient didjr^- 
These derivatives appear in (jS)), on account of (fTUj). In this respect, the equations are 
reminiscent of a system of quasi-linear diffusion equations. We shall discuss this matter 
further in Sectional 

ii) The frame variables E^^ enter into the remaining equations only through the spatial 
differential operator da- Recalling that is the spatial Hubble gradient we thus regard 
the variables Ea^ and as controlling the spatial inhomogeneity in the cosmological 
model, while the variables in Y directly determine the spacetime geometry (i.e. the 
gravitational field) and the matter content. 

2.4 The cosmological hierarchy 

We consider three classes of ever-expanding cosmological models, namely, Friedmann-Lemaitre 
(FL) models, spatially homogeneous (SH) models and general models, i.e. models without 
symmetry. The SH models admit a three-parameter local group G3 of isometrics acting on 
spacelike hypersurfaces, while in the FL models the isometry group is a Gg- For brevity we 
shall refer to the models without symmetry as Gq cosmologies. So we have the following 
cosmological hierarchy: 

FL C SH C Go. 

We will use S to denote the corresponding Hubble-normalized state spaces: 

iSpL C Ssu C Sgq- 

One can construct a more detailed hierarchy by including models with two or one spacelike 
Killing vectors. Classes of models with two Killing vectors are variously referred to'' as 
Gowdy^, T2-symmetric or G2 cosmologies. Models with one Killing vector are referred to as 
t/(l)-symmetric^ or Gi cosmologies. 

^See Uggla et al. 2003, Table 3, for selected references. 

^See Andersson et al. 2004 for a detailed analysis of Gowdy models using Hubble-normalized variables. 
^See Berger 2004 for a discussion of Mixmaster dynamics in these models. 
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SH cosmologies 

The SH cosmologies are obtained by requiring that the spatial frame derivatives of the 
gravitational field and matter variables Y, and of the normalization factor H be zero, i.e. 

doY = 0, r„ = 0. (15) 

It then follows that all the dimensional commutation functions and matter variables are 
constant on the hypersurface t = constant, which are thus the orbits of a three-parameter 
group 6*3 of isometries. The evolution equations (jH)) and Q imply that the SH restrictions 
(jl5|) define an invariant set of the full evolution equations, which we shall call the SH invariant 
set. Indeed, equation (jH)) is trivially satisfied, and equation (jH)) reduces to a system of ordinary 
differential equations, namely 

dtYA = Fa(Y). (16) 

The nontrivial constraints become purely algebraic restrictions on Y, which we write sym- 
bolically as 

C(Y) = 0. (17) 

An important consequence of this specialization is that the evolution equation (|7I) for 
decouples from the evolution equation for Y, which means that the dynamics of SH cosmolo- 
gies can be analyzed using only equations H^) and ff^ . In this context, one can think of 
the variable Y in equation © as defining a reduced Hubble-normalized state space, of finite 
dimension, for the SH cosmologies. The Bianchi classification of the isometry group then 
leads to a hierarchy of invariant subsets. 

In the SH context the restriction^^ f " = defines an invariant subset, giving the so-called 
nontilted SH cosmologies. This class of models have been analyzed in detail in the literature. 
We give an overview of the known results in Section 0] We note in passing that tilted SH 
cosmologies (i.e. f° 7^ 0) were first studied by King & Ellis 1973. Considerable progress has 
been made recently, using Hubble-normalized variables (see Hewitt et al. 2001, Hervik 2004, 
Coley & Hervik 2004). 

FL cosmologies 

The FL cosmologies are obtained by imposing the restrictions 

= 0, Va = 0, r„ = 0, 

-"^"This restriction means that the fluid 4-velocity is orthogonal to the orbits of the group of isometries 
(Ellis & MacCallum 1969). 
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which is equivalent to requiring that the fluid has zero shear, vorticity and acceleration^^. It 
follows that the density parameters are constant on the hypersurface t = constant 

In addition^^ the Weyl curvature is zero, 

and the 3-Ricci curvature satisfies 

5^/3 = 0, d^Qk = 0. (18) 

The dynamics of the FL models is governed by the matter evolution equations which 
simplify to a system of two ODEs, 

dtnm = [2q - (37 - 2)]n^ (19) 

dtnA = 2{q + l)nA, (20) 

with 

g = |(37-2)n„-fiA. (21) 
The Gauss constraint (fT^ simplifies to 

+ rik + = 1, (22) 

where all three quantities are functions of time only. For future reference we note that the 
curvature parameter Qk satisfies the evolution equation 

dtQk = 2gfifc, (23) 

as follows from (fT^ . ^ and (f^ . 

The remaining evolution equations simplify to 

dtEj = qEj, dtAa = qA^, dtNa/3 = qN^p, 

and hence decouple from ()19|) and (j^OJ. These variables reflect the choice of spatial frame, 
and do not affect the essential dynamics. Since q = q{t), these equations can be integrated 
to yield 

EJ = f{t)E^\ = f{t)A^, N^p = f{t)N^fs, 

where hatted variables depend only on the spatial coordinates x\ and f'{t) = qf{t). These 
quantities are restricted be equation (fTHj) . Models with negative, zero and positive spatial 

"'^"'^ Refer to WE, Section 2.4 for the characterizations of FL. 

""^^The first follows from the evolution equation for Va and the second from the definition of Ha and Tq. 
i^See Uggla et al. 2003, equations (AlO), (All), (A14) and (A15). 
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curvature are characterized respectively by 



fifc > 0, fifc = 0, fifc < 0, 

and in the corresponding canonical spatial frames we have 

Ar,^ = 0, A, = (A, 0,0), VLk = A\ 
No^p = 0, A^ = 0, fifc = 0, 
N^f, = diag(iV, iV, iV), A„ = 0, fifc = -\N\ 

Thus, the FL cosmologies with perfect fluid and cosmological constant are described by a 
reduced Hubble-normalized state space, namely the two-dimensional space with the state 
vector ^^a), and evolution equations (fT ^ - (PT|) . It follows from that for models 

with negative and zero spatial curvature {VLk > 0), the state space is bounded. This state 
space is shown in Figure [T] in Section |21 using VL^ and VL\ as variables. 

2.5 Cosmological equilibrium points 

In the analysis of the dynamics of SH models using the reduced Hubble-normalized state 
space iSsH, the equilibrium points (i.e. fixed points) of the dynamical system, defined by 

dtY = 0, 

naturally play a significant role. At this stage the Kasner vacuum solutions come into play, 
given in Hubble-normalized variables by 

iV,^ = 0, A, = 0, S2 = l, (]^ = 0, fiA = 0, g = 2, (24) 

with 

dtT.o.p = 0. 

These equilibrium points^^ form a 4-sphere (recall the definition (fT^ . and that S^/? is trace- 
free), which we shall call the Kasner sphere JC. If Sq,/3 is diagonal, then /C reduces to a 
circle. 

Two other important equilibrium points lie in the FL invariant set and satisfy 

Nal3 = 0, Aa = 0, = 0, Va = 0, 

^''The fact that the Kasner invariant set, defined by 124|) . consists only of equihbrium points is a conse- 
quence of the fact that we are using a Fermi-propagated spatial frame. See Uggla et al. 2003, page 7, for a 
different spatial gauge choice, which leads to a circle of equilibrium points, and "frame transition orbits" . 
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and one of tlie following sets of conditions: 

flat FL fi™ = l, ^^A = 0, g= 1(37 - 2) (25) 
de Sitter fi™ = 0, = 1, q = -1- (26) 

The local stability of the flat FL and de Sitter equilibrium points is of importance in con- 
nection with the phenomenon of isotropization, as discussed in Section 13.21 In the following 
tables we give the dimensions of the stable and unstable manifolds of these equilibrium 
points, for both nontilted and tilted SH models. 

Table 1: Stability of the flat FL and de Sitter equilibrium points in the state space for 
nontilted SH models (six dimensional) 

Equilibrium dimension of the dimension of the 

point stable manifold unstable manifold 

flat FL 2 4 

de Sitter 6 



Table 2: Stability of the flat FL and de Sitter equilibrium points in the state space for tilted 
SH models (nine dimensional) 

Equilibrium dimension of the dimension of the 
point stable manifold unstable manifold 



flat FL 5 4 

de Sitter, < 7 < | 9 

de Sitter, | < 7 < 2 6 3 

The SH equilibrium points, with the exception of de Sitter, represent self-similar solutions 
of the EFE, admitting a four-parameter group of similarities. We refer to WE, Section 
9.1, for a complete list of these self-similar solutions within the class of nontilted SH models. 

2.6 The silent boundary 

On account of ((Tj) and ((HI), the conditions 



Ej = 0, r„ = (27) 



^'^The de Sitter solution has other representations that involve the peculiar velocity See Tabled 
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define an invariant set, the so-called silent boundary^^, that forms part of the boundary of 
the Hubble-normalized state space. Solutions of the evolution equations in this invariant set 
do not correspond to solutions of the EFE, however, since the frame variables have to satisfy 
det(£'a*) 7^ 0. Nevertheless, it appears that this invariant set plays a fundamental role in 
describing the asymptotic behaviour of cosmological models. Indeed current investigations 
suggest that the past attractor and the future attractor for Gq cosmologies are subsets of 
the above invariant set (subject to A > in the case of the future attractor). We will 
discuss this matter further in Sections 121 and El The condition = in the definition 
of the silent boundary implies that d^Y = 0. It thus follows that the evolution equations 
and constraints, when specialized to the silent boundary, are precisely the SH evolution 
equations and constraints, as given by (fT^ and (fT7|). In addition each self-similar solution 
determines an equilibrium point in the silent boundary. 



3 Close-to-FL models and isotropization 
3.1 Close-to-FL models 

It is customary to regard a cosmological model as being "close to FL" if its metric is a 
linear perturbation of an FL metric (e.g. Bardeen 1980), i.e. the dynamics of such models are 
governed by the linearized EFE.^^ An alternative approach is to regard a cosmological model 
as being "close to FL" if an appropriate set of Hubble-normalized anisotropy parameters 
are small. For simplicity we will restrict our attention to two parameters. Firstly we 
consider the shear parameter E, defined by (fT^ . which describes the anisotropy in the rate 
of expansion of the fundamental congruence. Secondly we consider the Weyl parameter W 
defined by 

where Eap and Hai3 are the electric and magnetic parts of the Weyl curvature tensor relative 
to the fundamental congruence. The Weyl parameter can be thought of as quantifying the 
intrinsic anisotropy of the gravitational field. 

As described in Section 12.41 the FL cosmologies with negative or zero spatial curvature 
and positive cosmological constant have a two-dimensional Hubble-normalized state space. 
Using Qk and as independent variables, the state space is shown in Figure (H The 
equilibrium points represent the fiat FL model (F), de Sitter model (dS), and the Milne 
model (M). One sees that F is the past attractor and dS is the future attractor. The Milne 

^^This notion was first introduced in van Elst et al. 2002 and Uggla et al. 2003. See Andersson et al. 2004 
for further discussion. 

^^For a fully covariant approach to scalar perturbation of FL, which physically describe density fluctuations, 
we refer to Ellis & Bruni 1989, and Ellis et al. 1989. 

^^See WE pages 63-4, for a tentative definition of "close-to-FL" . Kristian & Sachs 1966 give a list of 
parameters, called "adjustable parameters" in their Table 1 (page 393) and referred to as "anisotropy and 
inhomogeneity parameters" on page 398. 
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point M is a saddle, which means that a subset of open FL models with A > will have a 
"close-to-Milne" epoch of finite duration. 




t V M 

open FL, A = 

Figure 1: The Hubble-normalized state space for flat and open FL cosmologies with A > 
and equation of state p = (7 — l)p, | < 7 < 2. 

When one regards the FL models as a subset of the SH models, a complication arises: a 
specific FL model is represented by infinitely many orbits in the reduced Hubble-normalized 
SH state space. This multiple representation occurs because a fiat FL model admits a G3 
of Bianchi type I and one of Bianchi type VHq, while an open FL model admits a G3 of 
Bianchi type V and one of Bianchi type Vllh, for any h > (see Ellis & MacCallum 1969). 
The state space can be described by the variables 

(^A, Qk, N), 

where 

nk = A'^ = hN^. 

Here h is the group parameter and A and are spatial curvature variables defined by 

N^p = diag(0, N,N), A, = (A, 0, 0) . 

The state space, which we will denote by 5fl, is shown in Figure |21 

We see that the state space is the union of a one-parameter family of "leaves" , and that 
the Milne model is represented by a line of equilibrium points Mh, one in each leaf. The 
family of orbits in one leaf with h > are physically equivalent to those in another. The 
orbits in the leaf with h = 0, and hence = 0, are multiple representations of fiat FL models 
with = or i^A > 0. In particular, the orbit with ^a = 0, = is an unbounded orbit 
that describes the fiat FL model, and the orbit with = 0, fim = is an unbounded orbit 
that describes the de Sitter model. The state space = is shown in Figure El 
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Figure 2: The Hubble-normalized state space for flat and open FL cosmologies with A > 
and equation of state p = (7 — | < 7 < 2, showing multiple representations. 

The FL state space <Sfl is an invariant subset of the SH state space Ssh, defined by the 
requirement that the shear of the congruence normal to the orbits is zero. It is natural to 
regard an SH model as being close to FL during some epoch if its orbit lies in a sufficiently 
small neighbourhood of the FL state space iSpL- The Hubble-normalized Weyl tensor is 
unbounded in such a neighbourhood, however, with the result that close-to-FL models can 
have large Hubble-normalized Weyl curvature. In particular there are models with 1, 
E and W large (close to flat FL) and models with JIa fa 1, E and W large (close to 
de Sitter), during some epoch. 

It is also of interest to consider the connection with perturbed FL models. Are orbits 
in a neighbourhood of 5fl described accurately by solutions of the linearized EFE, or do 
nonlinear effects come into play? Consideration of SH cosmologies of Bianchi type VHq 
shows that the latter occurs. Linearizing the Hubble-normalized evolution equations about 
the FL invariant set = and Qk — gives 

W ^ Ce=^(^-^)^ for I < 7 < 2, 

while an analysis of the exact equations for large r gives 

r Ce3(^-i)^ for I < 7 < I 
W ^ < CT-h\ for 7 = I 

[Cei(2-7)r^ for|<7<2. 
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de Sitter 



flat FL, A > 



— N 



flat FL, A = 



Figure 3: The state space — 0, showing multiple representations of the flat FL models. 

(Wainwright et al. 1999). Thus for this class of models, hnearization does not give the correct 
dynamics. 

It would be of interest to extend this analysis to Gq cosmologies, i.e. consider 5fl as an 
invariant set of the state space of Gq cosmologies and analyze the dynamics of close-to-FL 
cosmologies in this context. It would be necessary to determine all representations of the 
FL models within the Go state space, which would give an invariant set larger than »Sfl- 

3.2 Isotropization 

The term "isotropization" is used in three different contexts in cosmology: 

i) asymptotic isotropization into the future, 

ii) isotropic initial singularity, 

iii) intermediate isotropization. 

We discuss these dynamic phenomena in turn. 

Asymptotic isotropization into the future 

A well known result of Wald 1983 states that any SH cosmology with zero or negative spatial 
curvature (i.e. Qk > 0) approaches the de Sitter solution in the sense that 



Wald proved this result without making specific assumptions about the matter/energy dis- 
tribution. More detailed analyses of SH cosmologies with a tilted perfect fluid (Goliath & 



t— >+oo t— >+oo 



lim E = 0, lim Qa — 1. 
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Ellis 1999) revealed that a bifurcation occurs at 7 = |: if 7 > | the peculiar velocity 
of the fluid does not tend to zero. Specifically, lim Va = Ca, where c^c" < 1 if 7 = | and 

CaC°' = 1 if I < 7 < 2. This behaviour occurs for the following reason. For the de Sitter 
solution the reduced state vector Y in ® is given by 

N^p = 0, Aa = 0, S„^ = 0, Qm = 0, Qa = 1, (28) 

with Va subject to the evolution equation 

In the reduced Hubble-normalized state space the de Sitter solution is thus represented by a 
set of equilibrium points, given by (j2Hl) and with Va subject to Va = 0, or = 1 (a point and 
a sphere) if 7 7^ |, and subject to f ^ < 1 (a solid sphere) if 7 = |. The limiting behaviour of 
determines the future attractor in the reduced Hubble-normalized state space, as shown 
in Table H 



Table 3: The future attractor for SH cosmologies with A > and equation of state 
P = (7 - 1)P- 



Equation of state Future attractor 



0<7<| 


de Sitter with zero tilt {va -■ 


= 0) 


7 = 1 


de Sitter with intermediate 


tilt (0 < VaV'^ < 1) 


|<7<2 


de Sitter with extreme tilt 


(VaV'' = 1) 



The future attractor for SH cosmologies with A > forms the basis for the description of 
Go cosmologies that are future asymptotic to the de Sitter solution, as we now explain. The 
frame variables E^^ satisfy equation ((7|), which, when specialized to the de Sitter equilibrium 
solution (j28p can be solved to give 

jz/Qi -C/q e , 
where the E^^ are constants. It follows that 

lim Ej = 0. (29) 

In addition the constraints imply that the spatial Hubble gradient is zero (r^, = 0). Thus in 
the Hubble-normalized Gq state space, the de Sitter solution is described by orbits that are 
future asymptotic to a de Sitter equilibrium point on the silent boundary. 
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Go cosmologies that are future asymptotic to the de Sitter solution are described within 
the Hubble-normalized state space as follows. Along a fixed timeline, the evolution is de- 
scribed by an orbit that is future asymptotic to a de Sitter equilibrium point on the silent 
boundary. The specific equilibrium point will depend on the equation of state parameter 
7 according to Table El The idea is that the evolution is described by an orbit that is at- 
tracted to the silent boundary, i.e. holds, and within the silent boundary, the evolution 
is governed by the SH evolution equations, as described in Section 12.61 The asymptotic 
dependence of the Hubble-normalized variables for this class of models has been determined 
by systematically integrating the evolution equations (Lim et al. 2004). The fact that the 
orbit is asymptotic to the silent boundary means that the spatial dependence is asymptot- 
ically unrestricted, which manifests itself in the appearance of an arbitrary 3-metric in the 
asymptotic expansion of the spacetime metric (Lim et al. 2004). 

The issue of how large is the class of Gq cosmologies that are future asymptotic to de 
Sitter has not been completely resolved. For recent progress we refer to Rendall 2003. 

Isotropic singularities 

The second type of asymptotic isotropization, namely isotropization at the initial singularity, 
arises in connection with the notion of quiescent cosmology (cf., e.g.. Barrow 1978), which 
provides an alternative to cosmic inflation. The idea is that, due to entropy considerations 
on a cosmological scale, a suitable initial condition for the universe is that the Weyl curva- 
ture should be zero (or at least dynamically unimportant) at the initial singularity; this is 
the Weyl curvature hypothesis, proposed by Penrose 1979 (see page 630). This hypothesis 
leads to the notion of an isotropic initial singularity. Recently, Khalatnikov et al. 2003 (see 
page 3) have suggested that a solution with an isotropic singularity could represent an in- 
termediate asymptotic regime during expansion of the universe after inflation. Cosmological 
initial singularities of this type first arose as a special case in the general analysis of initial 
singularities performed by Lifshitz & Khalatnikov 1963 (see page 203). This type of initial 
singularity was also encountered in the work of Eardley et al. 1972 (see page 101). Subse- 
quently, Goode & Wainwright 1985 gave a formal definition of an isotropic initial singularity, 
using a conformally related metric, and derived various properties. Heuristically, the model 
approaches the fiat FL model near the singularity: 

lim S = 0, lim VLm = 1. 

t—>—oo f— >— oo 

For SH models, the occurrence of an isotropic singularity is characterized by the orbit in 
the reduced Hubble-normalized state space being past asymptotic to the fiat FL equilibrium 
point F. In other words, the orbit of a SH model with an isotropic singularity lies in the 
unstable manifold of the fiat FL equilibrium point, which confirms that isotropic singularities 
occur only for a set of initial conditions of measure zero. 

The above behaviour forms the basis for the description of Gq cosmologies that admit 
an isotropic singularity, as we now explain. The frame variables Ea satisfy equation (|7j). 
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which, when speciahzed to the flat FL equihbrium point can be solved to give 

E ' = E ^e^(^^~2)* 
where the E^ are constants. It follows that 

lim Ea' = 0, 

assuming | < 7 < 2. In addition, the constraints imply that the spatial Hubble gradient is 
zero {va = 0). Thus in the Hubble-normalized Go state space, the flat FL model is described 
by an orbit that is past asymptotic to an equilibrium point on the silent boundary. 

Go cosmologies that have an isotropic singularity are described as follows. Along a 
fixed timeline, the evolution is described by an orbit that is past asymptotic to the flat FL 
equilibrium point on the silent boundary. The asymptotic time dependence for this class of 
models has been determined by systematically integrating the evolution equations (Lim et 
al. 2004). The fact that the orbit is past asymptotic to the silent boundary means that the 
spatial dependence is asymptotically unrestricted, which manifests itself in the appearance 
of an arbitrary 3-metric in the asymptotic expansion of the spacetime metric (see Lim et 
al. 2004). 

Intermediate isotropization 

The flat FL equilibrium point F, given by (j^ . is a saddle point in the reduced Hubble- 
normalized state space for SH cosmologies. Thus, for any e > there is a family of orbits, 
corresponding to an open set of initial conditions, that pass through an e-neighbourhood 
of F but are not asymptotic to F. For the corresponding cosmological models there will be 
a finite time interval during which the rate of expansion of the model is highly isotropic 
(i.e. S ^ 1) and hence compatible with observational constraints. This behaviour is called 
intermediate isotropization., and occurs in models of all Bianchi types except Bianchi type I, 
in which case F is a sink.^^ 

At this stage, intermediate isotropization for Go cosmologies has yet to be investigated 
in detail. 

4 The dynamics of nontilted SH cosmologies 

The largest class of cosmologies for which there is a comprehensive and quite detailed knowl- 
edge of the dynamics, is the class of nontilted SH cosmologies with A > and equation of 
state p = (7 — 1)/), with | < 7 < 2. The reduced Hubble-normalized state space iSgH is 
an unbounded subset of M®. This class of cosmologies has two subsets referred to as class 
A and class B in the classification of Ellis & MacCallum 1969. We focus on the class A 
models, since they display a wider range of dynamical phenomena, while at the same time 

^^Bianchi type I cosmologies undergo asymptotic isotropization as t ^ +00. 
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permitting a simpler choice of frame. In describing the dynamics of this class of cosmologies, 
it is helpful to make use of the hierarchy of invariant subsets of lower dimension that arise 
from specializing the matter/energy content and the Bianchi type^°. We list these invariant 
subsets in Table El 

Table 4: Dimension of invariant subsets 



Invariant subset dimension 
zero cosmological constant I^a = 5 
vacuum flm = 5 
Bianchi VIq and VIIq 5 
Bianchi II 4 
Bianchi I 3 



We will discuss various dynamical phenomena, indicating how likely it is that each will 
occur. We say that a dynamical behaviour is 

i) generic, if it occurs for all initial conditions except for a set of measure zero^^, 

ii) typical, if it occurs for a set of initial conditions of positive measure, whose complement 
is also of positive measure. 

iii) special, it it occurs for a set of initial conditions of measure zero. 

One can say that generic behaviour has probability one, typical behavior has probability 
between zero and one, and special behaviour has probability zero. 

We now list the various dynamical phenomena in Table El and discuss each in turn. 

Table 5: Dynamical phenomena in nontilted SH cosmologies with A > and equation of 
state p = (7 — l)p, I < 7 < 2. 



Dynamical behaviour probability 
oscillatory singularity/Mixmaster dynamics generic 
inflationary isotropization/future asymptotic to de Sitter generic 
Weyl curvature domination typical 
intermediate isotropization typical 
isotropic singularity special 



^°The Bianchi types that occur in class A are I, II, VIo, VIIq, VIII and IX, with the last two being of 
maximum generality. We exclude Bianchi IX because models of this type can recoUapse. 

^^One can equivalently say "for all orbits in the reduced Hubble-normalized state space except for a set 
of measure zero" . 
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4.1 Oscillatory singularity /Mixmaster dynamics 

Within SsH, the Mixmaster dynamics is described by a past attr actor, which is the union 
of the Kasner circle of equihbrium points and the vacuum Bianchi II orbits, first described 
by Ma & Wainwright 1992^^. The attractor contains infinite heteroclinic sequences, i.e. infi- 
nite sequences of Kasner points joined by vacuum Bianchi II orbits. A generic orbit that is 
asymptotic to the attractor shadows one of these heteroclinic sequences. The evolution of the 
corresponding cosmological model towards the initial singularity is thus a non-terminating 
sequence of quasi-equilibrium Kasner states, punctuated by increasingly brief curvature tran- 
sitions. The early work on the attractor was based solely on heuristic local stability argu- 
ments and numerical simulations. Recently, however, Ringstrom 2001 proved^^ the existence 
of the past attractor, using the Hubble-normalized evolution equations. This result may have 
wider significance in view of the conjectures about the local past attractor for Go cosmologies 
that have recently been made (Uggla et al. 2003 and Section E)). 

As indicated in Tabled Mixmaster dynamics are generic. Examples of special behaviour 
are provided by models which are past asymptotic to a specific Kasner solution and by 
models with an isotropic singularity. 

4.2 Inflationary isotropization 

Within iSsH, inflationary isotropization is simply described by the fact that the de Sitter 
equilibrium point (|26|) . which describes the de Sitter solution, is a global sink^'', and hence 
forms the future attractor in the state space. This behaviour is generic, and special behaviour 
is restricted to the invariant set = 0. The dimensionless scalars satisfy 

lim (S, W, VLra) = 0, lim Qf, = 1. 
The asymptotic time dependence can be obtained from Lim et al. 2004 and Rendall 2003. 

4.3 Weyl curvature domination 

Generic orbits in the invariant subset = escape to infinity, and satisfy 

lim W = oo. (30) 

This result holds for Bianchi VIII orbits if the equation of state satisfies 1 < 7 < 2 (Horwood 
& Wainwright 2004) and for Bianchi VIIq orbits of 1 < 7 < 2 (Wainwright et al. 1999). In 
other words the Weyl curvature is dynamically dominant at late times. From a mathematical 
perspective this result means that the invariant subset Q/^ = does not admit a future 

^^This article is reprinted in Hobill et al. 1994. See also WE, Section 6.4. 

^■^Ringstrom's proof was given for the case of Bianchi IX cosmologies. A technical difficulty remains in 
extending the proof to Bianchi VIII cosmologies. 

^''This statement is essentially equivalent to the theorem of Wald 1983, discussed in Section 
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attractor. Despite this fact, one can still determine the asymptotic dependence of the state 
vector as t — > +00 (Horwood & Wainwright 2004). 

The asymptotic behaviour in the invariant set = as t — +00, as epitomized by (j30p . 
has a significant effect on the dynamics in the full state space fl\ > 0, in that it approximates 
the intermediate dynamics of a typical class of models with > 0. One can think of the 
asymptotic state of models with i^A = acting as a generalized saddle in the full state space, 
temporarily attracting orbits but eventually repelling them. In summary, for a typical class 
of models with I^a > there will be a finite epoch during which Weyl dominance occurs 
(i.e. >V> 1). 



4.4 Intermediate isotropization with ~ 1 

A cosmological model undergoes intermediate isotropization with Q^n ~ 1 if there is a finite 
epoch during which the model is close to the flat FL model. This behaviour is created by 
the local stability properties of the fiat FL model. We have seen that the flat FL model has 
a multiple representation, as the fiat FL equilibrium point F and as an unbounded orbit (see 
Section Em in particular. Figure E}- Orbits that are future asymptotic to F (i.e. orbits that 
form the two dimensional stable manifold of F) satisfy 

lim Qrn = 1, lim S = 0, lim W = 0. (31) 

t^+oo t~*+oo t—*+oo 

On the other hand orbits that are future asymptotic to the unbounded orbit satisfy^^ 

^ ^ ~ , (32) 
+00 if 1 < 7 < |. 

We say that models that satisfy ()31|) are asymptotic to flat FL in the strong sense, while those 
satisfy ()32|) are asymptotic to flat FL in the weak sense. Both of these asymptotic states 
can act as (generalized) saddles in the full state space, leading to intermediate isotropization 
with ~ 1- Thus intermediate isotropization with 0,^ ~ 1 in the full state space is a 
result of the occurrence of asymptotic isotropization with lim Q^n = 1 in lower dimensional 

t—*+ao 

invariant sets. 

Intermediate isotropization in the weak sense with Q^n ~ 1 is illustrated schematically in 
Figure |3 An orbit that is past asymptotic to a Kasner equilibrium point /C spirals around 
the fiat FL orbit, creating an epoch of intermediate isotropization with Q^n ~ 1, before 
approaching the de Sitter equilibrium point. 



4.5 Isotropic singularity 

Within SsYL, a cosmological model with an isotropic singularity is simply described by an 
orbit that is past asymptotic to the fiat FL equilibrium point (j^ . i.e. an orbit that lies in 

^^These orbits form the four dimensional Bianchi VIIq invariant set with J7a = 0. 
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Figure 4: Intermediate isotropization in the weak sense. 



the unstable manifold of this equilibrium point. Recall that this manifold is of dimension 4 
(see Table It follows that 



for these models. The asymptotic time dependence can be obtained from Lim et al. 2004. 

5 The vacuum past attractor and numerical simula- 
tions 

5.1 The vacuum past attractor 

We now describe the attractor in the Hubble-normalized Gq state space, following the dis- 
cussion in Uggla et al. 2003, but using a different choice of gauge. One of the goals of that 
paper was to give a precise statement of the so-called BKL conjecture, by describing the 
past attractor in the Hubble-normalized state space. 

The BKL conjecture. For almost all cosmological solutions of Einstein's 
field equations, a spacelike initial singularity is vacuum- dominated, local 
and oscillatory. 

For cosmological models with a perfect fluid matter source, the phrase "vacuum-dominated," 
or, equivalently, "matter is not dynamically significant," is taken to mean that the Hubble- 
normalized matter density Vtm tends to zero at the initial singularity. One might then expect 




lim (S, W, ^a) = 0, lim VL, 



'm 



1 
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that the past attractor for cosmological models with a perfect fluid matter source would be 
the same as the past attractor for an idealized vacuum cosmological model. It turns out 
(Uggla et al. 2003) that this expectation is unwarranted - the peculiar velocity plays a 
role in determining the past attractor. Nevertheless, partly in the interests of simplicity and 
also because the numerical simulations that have been recently performed (Garfinkle 2003) 
have been for vacuum models, we shall restrict our discussion here to the vacuum case. 

In terms of Hubble-normalized variables and the separable volume gauge, the spacelike 
initial singularity in a Gq cosmology is approached as t — > — oo. We now define a silent 
initial singularity to be one which satisfies 

lim Ej = 0, lim r„ = 0, (33) 

t— >— oo t^—oo 

lim daY = 0. (34) 

t—f — OO 

More precisely, we require that (jH^ and (jS^ are satisfied along typical timelines of Bq. One 
might initially think that the condition (jMj) is a consequence of ()33|) . since 

dY 

d Y = E ' — 

However the analysis of Gowdy solutions with so-called spikes^^ shows that the partial deriva- 
tives dY/dx"^ can diverge as t ^ — oo. Thus the requirement (jMjl demands that E^ tend 
to zero sufficiently fast. The asymptotic analysis of the Gowdy solutions shows that is 
satisfied along typical timelines even when spikes occur (see Uggla et al. 2003, Section IV. A). 
In more general G2 cosmologies, however, the validity of ()34|) is still open and further work 
is needed. 

We now refer to Uggla et al. 2003, pages 9-10, for further evidence to justify making the 
following conjecture. 

Conjecture 1: For almost all cosmological solutions of Einstein's field 
equations, a spacelike initial singularity is silent. 

Proving this conjecture entails establishing the limits (|33| ) - (j34j) . 

We think of the evolution of the Hubble-normalized state vector X()f:,x*) for fixed x*, 
as being described by an orbit in a finite-dimensional Hubble-normalized state space. As 
t ^ —00, this orbit will be asymptotic to a local past attractor A~, which in accordance with 
the definition of silent initial singularity, will be contained in the silent boundary, defined by 
(Elj). 

The past attractor is based on the Kasner solutions. We have seen that in the reduced 
Hubble-normalized state space for SH models, the Kasner solutions are represented by a 
sphere of equilibrium points. In the full Go state space, these solutions are represented by 
orbits that are past asymptotic to Kasner equilibrium points on the silent boundary. Thus 

26See Berger & Moncrief 1993, Garfinkle & Weaver 2003, Kendall & Weaver 2001. 
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we focus our attention on the Kasner sphere JC on the silent boundary, defined by 



Aa = 0, N^fs = 0, = 1, 

in addition to ()27p . Eacli non-exceptionaP^ equilibrium point on the Kasner sphere has one 
negative eigenvalue and hence has a one-dimensional unstable manifold into the past, which 
represents a SH vacuum solution of Bianchi type II on the silent boundary. We refer to these 
orbits as curvature transitions, since on each such orbit one degree of freedom of the spatial 
curvature is activated. These curvature transitions satisfy 

det(iV„^) = 0, := N^fsN'^^ - (iV,")' = 0, iV," ^ 0, 

corresponding to the fact that two eigenvalues of A^q,^ are zero. In addition 

Aa = 0, 

and the Gauss constraint simplifies to 

S2 + ^(iV„")2 = l. 

We now make our second conjecture, motivated by our knowledge of SH dynamics. 

Conjecture 2: The local past attractor for vacuum Go cosmologies 
with a silent initial singularity is 

A' = }CuTn, 

where /C is the Kasner sphere and T/v is the set of all curvature transitions 
in the silent boundary. 

Establishing this conjecture entails proving the following limits: 

lim = 0, (35) 

t^—oo 

lim det{Nap) = 0, (36) 

t—>—oo 

lim [N^pN^^ - (7V„")2] = 0. (37) 

t— »— oo 

The fact that the sequence of transitions is non-terminating implies that 

lim and lim N^'' 

t—>—oo t— >— oo 

^^The exceptional Kasner points describe the flat Kasner solution (the Taub form of Minkowski spacetime) 
and are characterized by the restriction E^^'^'E^^^ — Eq,^ = 0, where ( ) denotes traccfrce symmetrization. 

^^If one performs a spatial rotation to a shear eigenframe, one sees immediately from the evolution equa- 
tions for Nap that one diagonal component of N^p is unstable into the past. 
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do not exist. However it follows from the Gauss constraint ()12|) and equations ()33|). ()34|) and 
dnZj) that 

hm p2 + ^(iV,")2] = l. 

t— >— oo 

In addition it follows from one of the other constraints in (jllll . the Codacci constraints^, and 
the limits dM)), (El and (|n7|) that 

lim e"^^A/'«5S^^ = 0. (38) 

t—*—oo 

We thus expect that in a shear eigenframe, Na/3 will be "asymptotically diagonal" as t ^ —00. 

The geometrical interpretation of the asymptotic evolution is as follows: the orbit de- 
scribing the evolution of the state vector X(t, x*) for fixed approaches the Kasner sphere JC 
in the silent boundary and then shadows increasingly closely an infinite sequence of curvature 
transitions joining Kasner equilibrium points. 



5.2 Numerical simulations 

Some numerical simulations of the past asymptotic behaviour of vacuum Go cosmologies 
using the Hubble-normalized evolution equations (ffjl- (fTT|l have been recently performed by 
Garfinkle 2003. He found it convenient to use the deceleration parameter as a dynamical 
variable, with the defining equation (fTUj) forming an additional constraint. The evolution 
equation for q has the form of a quasilinear diffusion equation, and it is the only evolution 
equation to contain second order spatial derivatives. 

The Hubble-normalized evolution equations were solved numerically on a 3-torus, i.e. with 
periodic boundary conditions, using a Crank-Nicholson scheme. The initial conditions were 
chosen to be of the form 

EJ = ^lJ-^SJ, r„ = 0, iV„^ = 0, = -2i;-'Sjd,ij, 

(a2 cos y + as cos z + 62 + ^3 
ai cos x - as cos z &i - h 
—ai cos X — a2 cos y — &i — &2 

where arbitrary constants. 

The Gauss constraint gives a nonlinear elliptic PDE for ip{x'^): 

where is the Euclidean Laplacian in x, y and z. This PDE has to be solved numerically 
on a 3-torus to obtain explicit initial conditions. The remaining constraints are satisfied 
identically. Note that the above choice of corresponds to a conformally flat initial 3- 
metric. 

29See (2.28) in Lim et al. 2004. 




25 



The numerical simulations show a brief transient epoch after which a sequence of curva- 
ture transitions takes place, following the usual BKL transition law, which we now describe. 
The Kasner exponents Pa are related to the shear variables Tja/s in a shear eigenframe as 
follows: 

= diag(3pi - 1, 3p2 - 1, 3p3 - !)• 

The exponents satisfy 

and thus can be expressed in terms of a single parameter u according to 

u 1 + u u{l + u) 

The transition law between Kasner states can then be written in the form^° 

{u-1, ifu>2 
1 
, if 1 < M < 2. 
u — 1 

The numerical simulations also suggest that asymptotically Sq,/? and Najj have a common 
eigenframe in accordance with (|38|). 

The figures show the results of a numerical simulation with t = initially, and then be- 
coming negative, i.e. evolution towards the initial singularity {t —oo). Figures andlHt 
show the transitions of the Kasner parameter u along two different timelines and Fig- 
ures and inb show the diagonal values of A^^q,^ in the common asjTuptotic eigenframe, 
along the same timelines. Overall, the simulations provide support for the limits (j33|) and 
(jnH), and along typical timelines. 

The numerical simulation is incomplete in two ways. First, the duration of the Kasner 
epochs in the simulation is far too short, due to the fact that the decay of the eigenvalues 
of Nai3 to increasingly small values, which determines when the curvature transitions occur, 
is not described sufficiently accurately. The reason for this difficulty is that an arbitrarily 
chosen Fermi-propagated frame is not asymptotic to the eigenframe of Nap as t ^ — oo. 
Second, the simulation does not have high enough numerical resolution to correctly simulate 
spiky structures. These structures are known to develop in G2 cosmologies, on approach to 
the singularity (Berger & Moncrief 1993), and we expect that similar structure will develop 
in Go cosmologies. From our experience with numerical simulations of G2 cosmologies, these 
spiky structures occur within the particle horizon of timelines where one of the eigenvalues 
of Naf3 crosses zero (Lim 2004). As the particle horizon shrinks into the past, increasingly 
high numerical resolution is needed to simulate the spiky structures. It may be prohibitively 
expensive to perform such simulations in Go cosmologies. 

^"See, for example, WE page 236. 
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Figure 5: Transitions of the Kasner parameter u and the diagonal values of N^jS in the 
common asymptotic eigenframe, as t —>■ — oo along a fixed timeline, i.e. for fixed spatial 
coordinates x*. 
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Figure 6: Transitions along a different fixed timeline. 
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6 A normalization for recollapsing models 



The principal deficiency of Hubble-normalized variables is that they break down in an ex- 
panding cosmological model that reaches a maximum state of expansion {H = 0), as is 
possible in an FL model with positive spatial curvature. This difficulty can be addressed, at 
least for FL models and SH models of Bianchi type IX, by using a modified normalization 
factor first proposed by Uggla (WE Section 8.5.2). 

In the interests of simplicity we introduce the modified normalization procedure for the 
class of FL models with positive spatial curvature. For these models the normalization factor 
is defined by 

where is the curvature scalar of the metric induced on the hypersurfaces t = constant. 
Since % > 0, this factor is positive throughout the evolution, including at the time of 
maximum expansion H = 0. Dimensionless variables are defined in the usual way, using D 
instead of H: 

It follows that 

Qrn + = 1. 

Introducing a time variable f according to 

dt _ 1 
df^D' 

one obtains the following evolution equations for H and 

2(1 + q)HnA, 
where 

g = |(37- 2-371^a). 
We note that q is related to the usual deceleration parameter q according to 

q = H\ 

whenever H 0. The state space is bounded, being defined by the inequalities 

-1<H<1, 0<nf,<l, 

and is shown in Figure [7| The sign of H determines whether the model is expanding {H > 0) 
or collapsing {H < 0). The sign of q determines whether the model is decelerating (g > 0) 



dH 
df 
dflA 
~d^ 
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or accelerating. The fixed points of tlie evolution equations follows: 

CIa = 0, = 1, the fiat FL model F±, 

f^A = 1, = 1, the de Sitter model dS±, 

37-2 ~ 

f^A = , H = 0, the Einstein static model E. 

37 

The points F+ and dS+ represent expanding models (H = 1) while F_ and dS_ represent 
the time reversed models {H = —1). The orbits F_|_ — > F_ represent models that expand 
from a big-bang singularity and then recollapse to a future singularity. The orbits F+ dS+ 
represent models that expand indefinitely from a big-bang singularity, enter an accelerating 
epoch, and approach de Sitter at late times. 




Figure 7: State space for the FL cosmologies with positive and zero spatial curvature, A > 
and equation of state p = (7 — l)p, | < 7 < 2. 

Figure 13 illustrates the central role played by the Einstein static model within the class of 
FL cosmologies with positive and zero spatial curvature. The fixed point E is a saddle point 
whose stable and unstable manifolds are one-dimensional. The unstable manifold describes 
a one-parameter family of cosmologies that are past asymptotic to the Einstein static model, 
and hence do not have an initial singularity. These cosmological models have recently been 
proposed by Ellis & Maartens 2004 as a possible description of the early universe, a so-called 
emergent universe. 
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We note that this normahzation procedure can be generahzed to the class of nontilted 
SH cosmologies of Bianchi type IX. The normalization factor is (WE, equation (8.6)) 



where n^/j = diag^rin, n22, n^s) . The resulting state space, which gives a unified description 
of the dynamics of nontilted SH cosmologies of Bianchi type IX, VIIq, II and I, has not been 
explored in detail. In addition, a variation of this normalization has been given by Heinzle 
et al. 2004, and has been used to give a detailed analysis of the nontilted SH cosmologies of 
Bianchi type IX that are locally rotationally symmetric. 

7 Conclusion 

The dynamical systems approach has two key attributes. Firstly, it provides a geometric 
framework for analyzing the space of cosmological models. Secondly, it provides a system 
of first order evolution equations, whose specific form depends on the choice of gauge. This 
system of equations is potentially useful for giving a rigorous analysis of the asymptotic 
dynamics of cosmological models, and also for performing numerical simulations. The utility 
of this approach has been amply demonstrated as regards the SH cosmologies, and there are 
indications that it may be equally useful as regards the Go cosmologies. 
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